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We review results from FNAL E791 concerning D°-D mixing and CP violation in D meson decays. We have 
searched for mixing in semileptonic D° — > K + £~u decays and in hadronic D° — > K + -k~ and D° — > K + tv~ tt + it~ 
decays. We have searched for CP violation in D° — > K + K~ /tv + tv~ - and D+ -> <j>TT+/K K + / K+ K~ -k + / n+ -k + 
decays. Finally, we have measured the difference in decay widths AF between the two mass-eigenstates of the 
D°-D system. This parameter affects the rate of D°-D mixing. We combine our results with those from other 
experiments to obtain confidence intervals incorporating all published experimental data. 



1. INTRODUCTION 

FNAL E79lQ is a hadroproduction experiment 
studying the weak decays of charm mesons and 
baryons. The charm particles were produced by 
impinging a 500 GcV/c tt~ beam on five thin 
target foils. The most upstream foil consisted 
of platinum; the other foils consisted of carbon. 
All foils were separated by about 15 mm such 
that D mesons decayed predominately in the air 
gaps between foils. The experiment took data 
from September, 1991 to January 1992, record- 
ing the world's largest sample of charm decays at 
that time. The number of reconstructed events is 
over 200 000. With this data sample the experi- 
ment has studied charm production 0,^], charm 
lifetimes ||f|], rare and forbidden D decays ||, 
D°-D ° mixing [[| , CP violation [QJ|] , and several 
other topics. Here we focus on the following: (a) 

searches for D°-D mixing in semileptonic D° — > 
K + £~D decays and in hadronic D° — > K + tt~ 
and D — > K + Tr~ir + Tr~ decays; (b) measurement 
of the doubly-Cabibbo-suppressed decay D + — > 
K + tt~it + ; (c) search for CP violation in neu- 
tral D° — > K + K~ /n + ir~ decays and in charged 
D+ -> (hn + 



/K K+ /K+K-TT+ 



de- 



1 The collaboration consists of: CBPF (Brazil), Tel Aviv, 
CINVESTAV (Mexico), Puebla (Mexico), U. C. Santa 
Cruz, Cincinnati, Fermilab, Illinois Institute of Technol- 
ogy, Kansas State, Massachusetts, Mississippi, Princeton, 
South Carolina, Stanford, Tufts, Wisconsin, and Yale. 



cays; (d) measurement of the width difference AT 
between the two mass-eigenstates of the D°-D 
system. These results are published in Refs. |JJ^ 
|). Throughout this paper, charge- conjugate 
modes are included unless otherwise noted. 

The experimental apparatus consisted of a sil- 
icon vertex detector followed by a two-magnet 
spectrometer, two segmented Cerenkov count- 
ers for hadron identification, an electromagnetic 
calorimeter for electron identification, and iron 
shielding followed by scintillator counters for 
muon identification. The silicon vertex detector 
consisted of 17 planes of silicon and was used to 
reconstruct decay vertices downstream from the 
interaction vertex. The spectrometer consisted 
of 35 planes of drift chambers and two propor- 
tional wire chambers. The two dipole magnets 
bent particles in the horizontal plane and had 
p T kicks of +212 GeV/c and +320 GeV/c. The 
Cerenkov counters contained gases with differ- 
ent indices of refraction; together they provided 
n/K/p discrimination over the momentum range 
6-60 GeV/c. More details about the detector can 
be found in Ref. §. 

Data was recorded using a loose transverse en- 
ergy trigger. After reconstruction, events with 
evidence of well-separated interaction and decay 
vertices were retained for further analysis. Some 
of the main criteria used to select charm decays 
are listed in Table [l]. The most important crite- 
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rion is that of SDZ, defined as the distance be- 
tween the interaction and decay vertices divided 
by the error in this quantity. Values used for this 
criterion ranged from 8 (for D° and decays) 
to 20 (for longer- lived D + decays). 



Table 1 

Main criteria used to select D decays. 



Selection criteria 



Typ. value 



SDZ = 

(^dec ~ *i„t)/v^dcc + <J'L 8-20 

p T (transverse to D direction) < 250 MeV/c 

m ^ n \ Z dcc ~ Z targot cdgcl/ CT soc > 5 

D impact parameter 

w/r/t int. vertex < 60 ^m 

Xtrack < ^ 

t = m D x (z doc - z int )/p < 5 ps 



2. SEARCH FOR D°-D MIXING 

E791 has searched for D°-D° mixing via 
semileptonic D° — > K + l~v decays and hadronic 
D° — > K + ir~ and D° — > if+Tr^Tr+Tr - decays. 
For these searches we require that the D° be pro- 
duced via D* + — > D°tt + decay, and thus the fla- 
vor of the D° (or Z) ) when created is identified 
by the charge of the associated pion. The flavor 
of the D° when it decays is identified by the final 
state. With this information we measure the ratio 
r mix = T(D° -+ D° -» f)/T(D° /). Each type 
of decay studied (semileptonic or hadronic) has an 
advantage and a disadvantage: the Klv decays 
cannot be fully reconstructed due to the missing 
neutrino, and thus the decay-time resolution is 
degraded. The Kir/K iririr decays are fully recon- 
structed and thus have good time resolution, but 
they contain "background" arising from doubly- 
Cabibbo-suppressed (DCS) amplitudes that pro- 
duce the same final state. The DCS amplitudes 
do not contribute to the semileptonic decays. 



2.1. Semileptonic D° — > K+£ v Decays 

Semileptonic D° — ► K + £~P candidates were se- 
lected by requiring that there be a two-track ver- 
tex with SDZ > 8. One track was required to 
pass kaon identification criteria in the Cerenkov 
counters, and the other track was required to 
pass either electron identification criteria in the 
calorimeter or muon identification criteria in the 
scintillator counters following the iron shielding. 
Tracks identified as muons were required to have 
p > 10 GeV/c to reduce background from de- 
cays in flight. We define a quantity M min = 
p T + \Jp\ + where p T is the transverse mo- 

mentum of the Kl system with respect to the D° 
direction-of-flight (obtained from the interaction 
and decay vertex positions), and M Ke is the in- 
variant mass of the Kl pair. To reduce back- 
grounds, M min is required to be in the range 1.6- 
2.1 GcV/c 2 , and M K( is required to be in the 
range 1.15-1.80 GeV/c 2 . The upper cut on M Kl 
removes background from D° — > K~tt + decays in 
which the pion is misidentified as a lepton. After 
these cuts, the candidate D° is paired with a ■k ± 
track originating from the interaction vertex (and 
having p > 2 GeV/c) to form a Z)* ± . 

Since there is an undetected neutrino in the fi- 
nal state, the candidate D° momentum cannot be 
measured directly. However, using the direction 
of the D , the measured K and £ momenta, and 
assuming the parent particle mass to be that of a 
D , one can solve for the neutrino momentum up 
to a two-fold ambiguity. We use the solution re- 
sulting in higher D° momentum, as Monte Carlo 
(MC) studies indicate that this provides a bet- 
ter estimate of the true momentum. From MC 
studies we determine that the r.m.s. difference 
between the calculated and the true D° momenta 
is about 15%. 

To search for a mixing signal, we divide the 
electron and muon samples into "right-sign" (RS) 
and "wrong-sign" (WS) decays. The former have 
the charge of the kaon being opposite to that of 
the pion from the D* , whereas the latter have 
the charge of the kaon being the same as that 
of the pion. A WS decay would indicate D°-D 
mixing. To determine the numbers of events in 
the four samples (e and fi, RS and WS), we per- 
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form an unbinned maximum likelihood fit using 

, + ) and decay- 



<\d° 



the Q value ( 
time t (m D „ x (z dcc — z int )/p) for each event. 
For D* + — > D°tt + decays, the Q distribution 
is sharply peaked at 5.8 MeV. We also include 
in the fits the Q and t distributions of back- 
ground. The results for the electron sample are: 
N RS = 1237 ± 45 and N ws 



= 4.4™. There- 



suits for the muon sample are: N RS = 1267 ± 44 
and N ws = l.Stu.l- There is no indication of a 
mixing signal. Combining results from both Kev 
and K[iv samples gives r mix = (0.11 ^'2°)% or 



r mix < 0.50% at 90% C.L. 



2.2. Hadronic D° -> K+tz' 
K + tv~tv + tv~ Decays 



and D° 



D a — ► K + ir~ candidates were selected from a 
sample of two-prong decay vertices, and D° — > 
K + ir~TT + ir~ candidates were selected from four- 
prong vertices and from three-prong vertices with 
an extra track added. For the final event selec- 
tion, we use a two-layer neural network with 12 
input variables (for Kit) or 7 input variables (for 
Kmnr). These variables include the p T of the 
D° candidate with respect to the ir~ beam direc- 
tion, the p T of the D° with respect to the D° di- 
rection (obtained from the interaction and decay 
vertex positions), SDZ, the decay vertex fit x 2 , 
the track fit x 2 s, the Cerenkov counters' response 
for the K, etc. This selection results in eight sep- 
arate data sets: D° and D , Kit and Kttiytt final 
states, RS and WS decays. We subsequently per- 
form a single unbinned maximum likelihood fit to 
all data sets, constructing the likelihood function 
from the kinetic energy Q, the decay-time i, and 
the reconstructed mass m of each event. Back- 
grounds were carefully modeled and included in 
the fit. For each WS data set, the fit included con- 
tributions from a possible DCS amplitude. This 
amplitude results in a different t dependence than 
that due to mixing, and this allows one to par- 
tially discriminate between the two sources. The 
full expression for the WS t distribution (at small 
t where we have acceptance) is: 



dN ws /dt 



e- rt x 



(1) 



\A 



lAnixl 2 * 2 + 2Re(A BCS A* mix )t 



DCS I 



An example of the time dependences of the three 
terms is plotted in Fig. [l| 

The results of the fitting procedure depend 
upon whether we allow interference between the 
DCS and mixing amplitudes, and also whether 
we allow CP violation in any of the coefficients 
in Eq. (|l|). In the most general case of allow- 
ing CP violation in all coefficients, we obtain 

r mix (Z>° -_£>°) = (0.181°$ ±0.17)% and 

r mix( D ° -> = (0.70 ±0.18)%. Allowing 
CP violation in only the interference term gives 
r mix = (0.39 ±o'la ± 0-16)% or r mix < 0.85% at 
90% C.L. 



200 




0.5 1 1.5 2 
decay time (psec) 
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Figure 1. An example of the time dependence 
of D° — > K + tt~ decays due to the DCS ampli- 
tude (dashed), the mixing amplitude (dotted), 
and the interference between the two (dashed- 
dotted). The sum of all three contributions is 
solid. 



If we assume no mixing contribution (consis- 
tent with Standard Model predictions at our level 
of sensitivity), we obtain a rate for DCS decays. 



The results are: 



DCS 



{K-k) 



(0.68 t£!S ±0.07)% 
and r-acsiKit-Kit) = (0.25^;^ ±0.03)%. These 
values are approximately tan 4 9 C x (phase space), 
as expected, and are consistent with our mea- 
surement of the DCS charged decay D + — > 
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K + ir~ir + ||. For this latter measurement the 
final event sample is shown in Fig. ^. There are 
substantial backgrounds arising from misidenti- 
fied charm decays such as D° — ► K~tt + it + and 
D+ -> K+K-TT+ (both Cabibbo-favored). We 
simultaneously fit for these backgrounds and a 
D + — ► i^ + 7r _ 7r + signal, finding 59 ± 13 candi- 
date events in the peak. This gives a measure- 
ment r BCS (Kmr) = (0.77 ± 0.17 ± 0.' 




K 4 7t" n* invariant mass (GeV/c 2 ) 

Figure 2. The K + tt~tt + invariant mass spectrum 
for events passing final selection criteria. A peak 
resulting from the DCS decay D + — > K + ir~Tr + 
(and also from the singly-Cabibbo-suppressed de- 
cay Df — ► K + -k~tt + ) is clearly visible. 



to be largest. Because the incoming beam is 7r , 
the production cross section for D and D~ (in 
our acceptance) is a few percent larger than that 
for D° and D + ; this production asymmetry must 
be corrected for in order to discern a CP asym- 
metry. To do this we define the ratios: 



Vd- 



N D ^ f /N D _ 



Then A 
if e, 



CP 



-K-tt+(tt+) 
K+n-(n-) 



(3) 
(4) 



n 



D^f/ £ D^K-ir+(Tv+) 



£— 



r) 



le- 

»/' D^K+ir-(ir~y 



where Sj^^j is the overall detection efficiency for 
D — > / (including acceptance). This relation- 
ship among efficiencies holds well for E791. We 
assume there is negligible CP violation in the 
Cabibbo-favored decay modes used to normalize 
the Cabibbo-suppressed decay rates [Eqs. (||) and 



3.1. Neutral D° Decays 



We measure A CP for D° 



K+K- and D° -> 
where the D° originates from D* + — > 
D°tt + and thus the flavor of the D° is identified 
by the charge of the associated pion. The final 
mass plots are shown in Fig. ^ along with those 
for the normalization channel D — > K~n + . The 
solid curves superimposed are our fits to the his- 
tograms. The integrals of the Gaussian distribu- 
tions used for the signals determine the number 
of D° — ► Ktt/KK/tttt decays. These event yields 
(listed in Fig. |J) are used to calculate t] d0 and 
7]— and subsequently determine A cp . The re- 
sults are: 



3. SEARCH FOR CP VIOLATION 

E791 has searched for CP violation in both 
charged and neutral D decays. For these searches 
we measure the time-integrated asymmetry A CP , 
defined as: 



A, 



CP 



T{D -> f) - TjD -> /) 
T(D - /) + T(D /) 



(2) 



We study only Cabibbo-suppressed final states, as 
for these modes CP-violating effects are expected 



A CP (K+K~) = 
A cp (tt+tt-) = 



-0.010 ± 0.049 ± 0.012 (5) 
-0.049 ± 0.078 ± 0.030. (6) 



These correspond to 90% confidence intervals 
-9.3% < A CP {K+K-) < 7.3% and -18.6% < 
A cp {it + tt-) < 8.8%. 

3.2. Charged D+ Decays 



We measure A CP for D + decays into the fi- 
nal states 4>tt+ ((f) -> K+K-), K*°K+ (K*° -► 
K~tt + ), K + K~tt + (nonresonant), and 7r + 7r~7r + . 
For all modes, the normalization channel is D + — > 



5 



800 
600 
400 
200 




D° -> K"7T* 
6075 ± 87 



_L 



1.7 1.8 1.9 2 2.1 2.2 

M(K~tt + ) (GeV/c 2 ) 




1000 
750 
500 
250 


60 
40 















D° - 


-> K*7T" 






719 


8 ± 95 











1.7 1.8 1.9 2 2.1 2.2 

M(K + tt") (GeV/c 2 ) 









1 D° -> K*K" 




ij 333 ± 21 







1.7 1.8 1.9 2 2.1 2.2 

M(K + K") (GeV/c 2 ) 



1.7 1.8 1.9 2 2.1 2.2 

M(K + K") (GeV/c 2 ) 




1.9 2 2.1 

M(tt + tt") (GeV/c 2 ) 



1.9 2 2.1 

M(tt + tt") (GeV/c 2 ) 



Figure 3. Final mass plots for D° (left) and 
D (right) decays into final states Ktt (top row), 
K + K~ (middle row), and 7r + 7r _ (bottom row). 



For the </>7r + final state we require 
\ m K+K- - m <f>\ < 6 MeV/c 2 ; for the K*°K+ final 
state we require \m K - v+ — m K , \ < 45 MeV/c 2 . 
The mass spectra for the final event samples are 
shown in Fig. EL and the measured asymmetries 
A CP are listed in Table ||. We also list measure- 
ments from other experiments, and for each decay 
mode we calculate a 90% confidence interval for 
A CP incorporating all measurements listed. As 
the measurements are from independent experi- 
ments, we assume their statistical and systematic 
errors uncorrelated. We observe that A CP for 
D + — > K + K~ir + is relatively well-constrained: 
the 90% CL interval is -1.6% to +2.0%. For 
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Figure 4. Final mass plots for D + (left) and D 

(right) decays into final states 4>tt, K K, KK-k 
(nonresonant), and tttttt. 
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Table 2 

CP asymmetries measured for D° and D + decays. Also listed are measurements from other experiments, 
and listed in boldface type are 90% confidence intervals incorporating all the experimental results. 



Mode 


A cp (E791) A cp (Others) 


D° -> K+K- 


+0.024 ± 0.084 (E687 @) 
-O.OfO ± 0.049 ± 0.012 +0.080 ± 0.061 (CLEO 0) 

-0.001 ± 0.022 ± 0.015 (E831 g§) 
-2.6 < A cp < 4.4 % 




-0.049 ± 0.078 ± 0.030 +0.048 ± 0.039 ± 0.025 (E831 g§) 
-4.1 < A cp < 9.2 % 


D+ -> K+K-TT+ 


-0.031 ± 0.068 (E687 M) 
-0.014 ± 0.029 ^ 

+0.006 ± 0.011 ± 0.005 (E831 @) 

-1.6 < A cp < 2.0 % 


D+ -> 4>TT+ 


-0.028 ± 0.036 0.066 ± 0.086 (E687 @) 
-6.8 < ^ cp < 4.1 % 


D+ -> K*°K+ 


-0.010 ± 0.050 -0.12 ± 0.13 (E687 @) 
-10 < A cp < 5.2 % 


D + — > 7r + 7r _ 7r + 


-0.017 ± 0.042 

-8.6 < A cp < 5.2 % 



4. MEASUREMENT OF THE WIDTH 
DIFFERENCE AT 

E791 has measured the difference in decay 
widths between the two mass-eigenstates of the 
D°-D system. This provides a measurement of 
the mixing parameter y = AT/(2T). Theoreti- 

cally, r mix = T(D° -> D° -> /)/r(£>° -_/) = 
(x 2 + y 2 )/2 for a;,?/ small, where a; = Am/T. 

The method used to determine Ar is as fol- 
lows: assuming no CP violation, the two mass- 
eigenstates are CP eigenstates and can be writ- 
ten D 1 = [\D°) + \D°))/y/2 and D 2 = (|D°) - 
\D°))/y/2. Observing a K + K~ final state 
(CP = +1) denotes a D 1 — > K + K~ decay, and 
dN KK /dt oc e -1 " 1 *. Observing a if~7r + or K + ir~ 

final state denotes a D° oi D° decay, respec- 
tively (neglecting DCS amplitudes for simplicity). 
In this case dN KlT /dt oc e" rt cosh(Ar/2)t (H), 
where T = (r x + T 2 )/2 and AT = r x - T 2 . Our 



previous limit r mix < 0.50% implies \y\ < 0.10 
or |Ar| < 0.48 ps"\ and thus cosh(Ar/2)£ » 1 
for the range of lifetime t accepted by the ex- 
periment. Thus dN KlT /dt cx e _rt , and Ar = 
2(T KK -r Kv ). Equivalently, y = AT/(2T) = 
t Kit /t kk — 1. From an experimental point of 
view, it is convenient that cosh(Ar/2)i w 1 as 
we actually measure t — t cut = t', and while 
e~ rt oc e~ rt , cosh(Ar/2)i' is not proportional 
to cosh(Ar/2)t. 

The D" -> K-TT+ and D° -> i^+Jf" final 
event samples are shown in Fig. |5|. The re- 
sultant lifetime distributions (for t') are shown 
in Fig. ^. These distributions have the back- 
ground shape subtracted, and fitting to them 
yields: T Kv = 2.420 ± 0.019 ps" 1 and T KK = 
2.441 ± 0.068 ps^ 1 . The difference in widths Ar 
is 0.04 ± 0.14 ± 0.05 ps _1 , where the first error is 
statistical (resulting from the fits) and the second 
error is the sum in quadrature of the systematic 
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errors. These systematic errors are listed in Ta- 
ble ||. The result corresponds to a 90% confidence 
interval -0.20 < Ar < 0.28 ps" 1 . For y we ob- 
tain 0.008 ±0.029 ± 0.010 or -0.042 < y < 0.058 
at 90% C.L. Combining this with a recent mea- 
surement by FNAL E831 (y = 0.0342 ± 0.0139 ± 
0.0074 HI) gives 0.006 < y < 0.052 at 90% C.L. 



12000 
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Figure 5. Final mass plots for D° — > K tt + (top) 
and L>° ->■ K + K~ (bottom). The peak to the 
right of the D° — > K + K~ peak (bottom plot) is 
due to D° — ► K ~ir + decays in which the pion has 
been misidentified as a kaon. 



5. SUMMARY 

We have searched for D°-D° mixing in both 
semileptonic and hadronic D° decays and see no 
evidence for it. We subsequently set 90% CL 
upper limits r mix < 0.50% and r mix < 0.85%, re- 
spectively. Assuming the mixing amplitude to 
be negligibly small (as predicted by the Stan- 
dard Model), we fit our data for the rate of DCS 
decays D° — > K + ir~ and D° — > K + 7r~7r + ?T~, 
obtaining r DCS {Kir) = (0.68 ±g;|| ±0.07)% and 

r Dcs(^ 7r7r7r ) = (°- 25 -o.u ±0.03)%. These re- 
sults are consistent with the rate we measure 
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Figure 6. Reduced lifetime (t 1 ) distributions for 
D° -> R-TT+ (top) and D° -> K+K~ (bot- 
tom). The contributions from background have 
been subtracted. 



for the DCS charged decay D + — > K + ir tt + : 
r BCS {KTTTr) = (0.77 ± 0.17 ± 0.08)%. 

We have also searched for CP violation in neu- 
tral D° — > K + K~ /ir + Tr~ decays and in charged 

D+ (t>TT+ /K*°K+ /K+K-TX+ /ir+TT-n+ de- 
cays. We see no evidence for CP violation, and for 
each mode we set 90% C.L. limits on the asym- 
metry parameter A CP . These results are listed in 
Table |. 

Finally, we measure the difference in decay 
widths Ar between the two mass-eigenstates of 
the D°-D system. We convert this into a mea- 
surement of the mixing parameter y = AT/(2T). 
Our result is y = 0.008 ± 0.029 ± 0.010 or 
-0.042 < y < 0.058 at 90% C.L. Combining this 
with a recent measurement by FNAL E831 gives 

0. 006 < y < 0.052 at 90% C.L. 
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